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A Method of Providing a Signal Extension for a Finite Length Signal and 

Apparatus Therefor 

Field of the Invention 
This invention relates to a method of providing a signal extension to a 
finite length signal and to an apparatus therefor, especially, though not 
exclusively to a method of providing a signal extension to a wavelet 
transformed signal of finite lengt, for example, in an image compression system. 
In particular, the invention can be utilised in an image compression system for 
performing subband coding of large images partitioned (segmented) into non- 
overlapped tiles. 

Background of the Invention 
One of the most efficient techniques for image compression is subband 
coding, such as wavelet transformation, where an image input signal is 
decomposed into several frequency subbands that are quantised and then coded 
for onward transmission or storage. Since this process is reversible, the original 
image signal can be reconstructed back from the coded subband information. 
An important practical problem in subband signal coding is that a high degree 
of complexity is introduced by decomposition /reconstruction when the signal 
length is substantial. This is especially critical for large images. The most 
commonly used approach to cope with this problem is to break an image into 
sub-images, known as blocks or tiles, of a smaller size and process them 
independently. This method, which is usually called non-overlapped tiling, can 
significantly reduce the complexity, but introduces very noticeable boundary 
artifacts in the regions of the edges of the tiles in the reconstructed image. Since 
a human eye is very sensitive to edges especially in uniform areas, these types of 
distortions are very noticeable and consequently degrade image quality 
dramatically. The problem of smooth boundary reconstruction is, of course, not 
limited to image coding, but exists in different aspects in many other 
applications. 

Due to its practical importance, the problem of reducing boundary 
artifacts has been extensively investigated and described in many publications. 
One known approach is to use overlapping samples from adjacent tiles, as 
suggested by II Kye Eom, Yoon Soo Kim and Jae Ho Kim in an article entitled 
"A Block Wavelet Transform for Sub-image Coding /Decoding" published in 
SPIE Vol. 2669 at pagesl69-177 and in US Patent No. 5,710,835. Another known 
approach is to employ a variety of post processing techniques at the image 
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reconstruction stage, as described by B. Jeon and J. Jeong in an article entitled 
"Blocking Artifacts Reduction in Image Compression with Block Boundary 
Discontinuity Criterion" published in IEEE Trans, on Circuts and Systems for 
Video Tech. Vol 8, N3, 1998 at pages 345-367 and by J. K. Su and R. M. 
5 Mersereau in ana article entitled "Post-processing for Artifact Reduction in 
JPEG-compressed Images" published in ICASSP-95, Vol. 4 at pages 2363-2366. 
Unfortunately, most existing methods can only be applied in practical systems 
with great difficulty due to the substantial additional complexity which they 
introduce. 

10 In this specification, including the claims, the terms "comprises", 

"comprising" or similar terms are intended to mean a non-exclusive inclusion, 
such that a method or apparatus that comprises a list of elements does not 
include those elements solely, but may well include other elements not listed. 

Brief Summary of the Invention 
15 The present invention therefore seeks to provide an image compression 

system and to a method of coding therefor which overcomes, or at least reduces 
the above-mentioned problems of the prior art. 

Accordingly, in a first aspect, the invention provides a method of 
providing an extension to at least one end of a signal of finite length at at least 
20 one boundary thereof, the extension being formed by reproducing a portion of 
the finite length signal adjacent to the boundary in an point symmetric fashion 
about a point on the finite length signal at or close to the boundary. 

In a second aspect, the invention provides a method of providing an 
extension to at least one end of a signal of finite length, the extension being 
25 formed by the steps of: 

defining a point at the end of the finite length signal; 

determining a length of the finite length signal starting from the defined 

point; 

duplicating the determined length in a point symmetric fashion about the 
30 defined point so as to provide an extension of the finite length signal beyond the 
end. 

In a preferred embodiment, signal extensions are provided at both ends 
of the finite length signal. It will be appreciated that the finite length signal 
could be image data, or could be other data such as speech or other acoustic 
35 data, or other data. 

According to a third aspect of the present invention, there is provided a 
method of subband image decomposition comprising the steps of: 
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receiving image data for an image that has been divided into a plurality 
of image blocks; 

applying a method of providing an extension to at least one end of a 
signal of finite length, as described above, to the image data for each image 
5 block; and 

decomposing the extended image data for each image block to provide 
subband decomposed image data. 

Preferably, the step of decomposing the extended image data comprises 
transforming the extended image data using a transform, such as a wavelet 
10 transform. In one embodiment, at least part of the decomposed image data is 
replaced by data formed from a combination of part of the image data and part 
of the decomposed image data. 

In a fourth aspect, the invention provides a method of compressing 
images comprising the steps of: 
15 receiving image data representing at least one block of an image which 

has been divided into blocks; 

decomposing the received image data using the method of subband 
image decomposition as described above; 

quantising the subband decomposed image data; and 
20 coding the quantised subband decomposed image data to provide 

compressed image data. 

In a fifth aspect, the invention provides a method of image tile 
reconstruction comprising the steps of: 

receiving at least two subband sequences of data; 
25 applying the method of providing an extension to at least one end of a 

signal of finite length at at least one boundary thereof as described above to the 
subband sequences of data; and 

performing subband synthesis on the extended subband sequences to 
produce reconstructed image tiles. 
30 In a preferred embodiment, two subband sequences are received. In one 

embodiment, at least one of the received subband sequences includes combined 
data formed from a combination of original image data and decomposed image 
data. Preferably, the step of performing subband synthesis includes the step of 
combining at least the received combined data with other data from at least one 
35 of the received subband sequences. 

Brief Description of the Drawings 
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One embodiment of the invention will now be more fully described, by 
way of example, with reference to the drawings, of which: 

FIG. 1 shows a simplified block diagram of a well-known two-channel 
one-dimensional analysis and synthesis filter bank; 
5 FIG. 2 shows a sequence of data samples in a sequence of finite length 

having extensions produced using a common symmetric method; 

FIG. 3 shows graph of an example of a signal across an image block, 
together with extensions generated using the conventional symmetric method 
and the new point-symmetric method according to the invention; 
10 FIG. 4 shows sequences of data samples using point-symmetric extension 

for whole-sample and half-sample types of symmetry; 

FIG. 5 shows a graph of a signal adjacent one of its boundaries to 
illustrate the differences between the point-symmetric and conventional 
symmetric extensions; 
15 FIG. 6 shows sequences of data samples to illustrate subband analysis 

and synthesis of odd-length sequences using odd-length filters; 

FIG. 7 shows sequences of data samples to illustrate subband analysis of 
even-length sequences using odd-length filters; 

FIG. 8 shows a block diagram of an apparatus for performing one- 
20 dimensional subband decomposition of sequences according to an embodiment 
of the present invention. 

FIG. 9 shows a block diagram of an apparatus for performing one- 
dimensional subband reconstruction of sequences according to an embodiment 

of the present invention. 
25 Detailed description of the invention 

Subband decomposition can be applied to images in a number of 

different ways depending on the application. Never the less, most of them are 

based on the same analysis /synthesis operation schematically presented in Fig. 

1. In general, there are two analysis filters in a first, analysis section: a low pass 

30 filter H 0 2 and a high pass filter Hj 3 and a synthesis section, which also includes 
two filters: a low pass filter G 0 6 and a high pass one G^7. An input sequence 
X[n] from an input 1 is low and high pass filtered firstly by filters 2 and 3. Then 
the low and high pass components from the filters are downsampled by a factor 
of two, producing outputs L[n] and H[n]. The low frequency subband L[n] and 

35 the high frequency subband H[n] can be used for further decomposition or 
quantised. At the synthesis stage, the low frequency subband L[n] passing 
through flow 4 and the high frequency subband H[n] passing through flow 5 
are upsampled by factor of two, low and high pass filtered by low pass filter G 0 
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6 and high pass filter G a 7. Finally, both subbands are summed by adder 8 to 
produce an output sequence Y[n] at an output 9, which would be equal to X[n] if 
the operation has the property of perfect reconstruction. 

A broad class of linear filters such as QMF, CQF and wavelets can be 
5 used in H 0/ H a , G 0 and G r However, regardless of the filter type, a proper 

extension of the input sequences X[n], L[n] and H[n] must be applied to perform 
the convolution properly. Very good results can be achieved by using the 
known symmetric extension shown in Fig.2. It is based on a mirror symmetry 
and the supplementary samples x[i] (dashed lines) can be easily found by 
10 copying the signal samples X[i] (firm lines) around an axis of symmetry 11 and 
13 as illustrated by arrows 15. This operation can be expressed by a formula: 
x[i] = X[i] , 

where i = [1, k] for odd-length filters, or i=[0, k-1] for even-length filters, 
where k is the number of samples required to extend the sequence. 
15 Even though this extension is more efficient than a known periodic one 

(where the signal is simply repeated at the extension) and is currently employed 
in most subband systems, it causes boundary artifacts when images are tiled, as 
shown in Fig. 3, where the narrow solid line shows the signal over a complete 
tile (Tile 2) and adjacent portions of adjacent tiles (Tiles 1 and 3). The bold line 
20 shows the reconstructed signal reconstructed using the known symmetric 
extension and shows how the reconstructed signal is discontinuous at the 
boundaries of the tiles. The dotted line shows the reconstructed signal 
reconstructed using the point symmetric method according to an embodiment 
of the invention and clearly shows that the discontinuity at the boundaries is 
(3 25 eliminated, or at least substantially reduced. This disadvantage of the 

conventional symmetric approach is also shown in Fig. 2, where there is no 
smooth transition between the real and extended sequences. Instead, sharp 
peaks may be produced, as around X[0] for example, that badly affect 
compression. 

30 To overcome the problems, an embodiment of the present invention 

proposes generating the extensions using a different type of symmetry, called 
hereinafter "point-symmetry" since it uses a different type of symmetry from the 
single-mirror-symmetry of the conventional extension. An example of the point- 
symmetric extension is shown in Fig. 4. As shown, the supplementary samples 

35 x[i] (dashed lines) are produced from the signal samples X[i] (firm lines) by 
arrows 27 going through a point of symmetry 25. Actually, two types of the 
point symmetry are possible: whole-sample (WS) and half-sample (HS). In the 
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first one, a point of symmetry 25 belongs to a boundary samples 17 or 19. In the 
second one, the point of symmetry 25 belongs to an imaginary sample 21 or 23 
located at the half way between others. Thus, the difference between those types 
of point symmetry can be considered as the difference in the point of symmetry 
5 25 positioning. 

In comparison with the conventional method the proposed one: 

- minimises the deviation at the tile boundaries of the reconstructed sequences, 
thus reduces boundary artifacts; 

- produces more smooth joint between the signal and extension by eliminating 
10 the slope discontinuity, thus leads to compression improvement. 

The efficiency of reconstruction with the new extension is illustrated in 
Fig. 3. As shown, the use of the point-symmetric extension results in more 
accurate reconstruction at the tile boundaries. Thus, it produces a smaller step 
between adjacent segments of the reconstructed signal approaching that of the 
1 5 originally smooth signal. 

A simple formula can be used to calculate the extension x[n] at the 
beginning and the end of X[n]: 

x[i] = 2*P-X[il, (2) 
where i = [1,..., k] for odd-length filters, or i=[0, k-1] for even-length filters 

20 P - value of the starting /terminating point of symmetry. 

If the beginning of X[n] is extended, the starting point of symmetry 
belongs X[0]. The terminating point of symmetry is the last sample of X[n] (X[7] 
in our example shown in Fig. 4 ). As for the number k of the extension samples 
required, it is equal to the number of the filter coefficients falling outside the 

25 sequences when the filter is positioned at the starting /terminating sample. 

Although this rule is similar to that one for conventional extension, the required 
number of the extension samples can be smaller, because the 
starting/terminating sample can be located not at the tile boundary. This will be 
illustrated by an example below. 

30 The difference between two extensions applied to the same signal is 

illustrated in Fig. 5. One may notice that the conventional method can produce 
an artificial sharp peak at the signal boundary, while the point-symmetric 
extension creates a very smooth joint by eliminating the slope discontinuity. 

The same as the conventional method the proposed approach operates 

35 slightly different for odd /even-length filters and odd /even-length sequences. 
Firs of all, WS is applicable to odd-length filters, while HS is used with even- 
length filters. An example of data flow in odd-length analysis/synthesis filter- 



W 

-7- 

bank applied to an odd-length sequence is presented in Fig. 6. At first stage the 
input sequence X[n] is point-symmetric extended (Fig. 6a). After applying the 
subband analysis, two sequences L[n] (Fig. 6b) and H[n] (Fig. 6c) are produced. 
The overall number of subband samples ( L[n] and H[n] ) is equal to the number 
5 of input samples and L[n] subband always contain one more sample than H[n], 
At the synthesis stage (Fig. 6 d-e) the extension samples are produced using the 
same rule, but in the high-band the point of symmetry is equal to 0 (Fig. 6e). The 
result of the analysis Y[n] (Fig. 6f) has the same number of samples as X[n]. It 
should be pointed out that this method for many wavelet filters produces L[0] = 

10 X[0] and L[M] = X[N], where M and N are maximum indexes of L[n] and X[n] 
correspondingly. This feature permits to reduce the number of convolution 
operations in the low-band channel by two and reduce the number of the 
extension samples required. 

In the case of even-length sequences, data reordering is required to 

15 preserve the perfect reconstruction property of the filter bank. Although 

reordering can be carried out in a number of ways it always aims to add missing 
information about one of the P to the output. An example of data flow in odd- 
length analysis/ synthesis filter-bank applied to an even-length sequence is 
presented in Fig. 7. Due to the point-symmetric extension of X[n], the last 

20 sample in the high-band is always equal to 0 (see Fig. 7c). Thus, its position can 
be used to store some additional information required for lossless 
reconstruction. One possible solution is illustrated in Fig 7d. H[3] is produced 
by subtraction L[3] from X[7] ( arrow 120). Then L[3] in the low-subband is 
replaced by X[7] (see Fig. 7e, arrow 122). After this rearrangement both 

25 subbands contain the same number of samples. At the synthesis section 
backward data reordering is required. At first, L[3] is set as the point of 
symmetry of the L[n] end (see Fig. 7g, arrow 124). Then L[3] is replaced with 
the sum L[3]+H[3] (see Fig. 7g, arrow 126) and H[3] set to 0 as shown in Fig. 7f. 
Finally, the sequences L[n] and H[n] are extended, filtered and summed 

30 together to produce Y[n]. If quantisation is not applied to L[n] and H[n] the 
output Y[n] = X[n], because this method preserves the perfect reconstruction 
property. 

The block diagrams of apparatus implementing the subband 
decomposition is shown in Fig.8. From the information about low-pass filter 
35 length 12 and high-pass filter length 14 applied to the first and second input of 
the extension samples counter 16, it calculates the number of extension samples 
k as described earlier. Then k coupled with the extension type (WS if the filters 
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are odd- length, else HS) are passed to the first input 18 of the extension samples 
generator 26. Two P values corresponding to the starting and terminating points 
of symmetry are passed from the point of symmetry selector 24 to the second 
input 20 of the extension sample generator 26. The selector 24 copies those 
5 values from the starting and the terminating samples of the input sequence X[n] 
from input 10, which is also sent to the third input 22 of the extension samples 
generator 26. The extension samples x[n] are calculated in the extension samples 
generator 26 using formula (2). Those samples are sent to the second input 32 of 
the analysis filter bank 30 and together with input samples from 10 applied to its 

10 first input 28 are used to perform the subband decomposition as illustrated in 

Fig.6b-c and Fig.7b-c. The output of the analysis filter bank, which consists of the 
low-pass band 34 and the high-pass band 36 is passed to the first and second 
inputs of the data reordering unit via the first and second output of the analysis 
filter bank correspondingly. If the signal length applied to the third input 38 of 

15 the data reordering unit 40 is even, the data reordering is performed as it has 
been described. One example of such reordering is shown in Fig. 7d-e. If the 
signal length is odd, the data reordering unit 40 simply passes by the samples 34 
and 36. The sequences L[n] and H[n] from the first 42 and second 44 output of 
the data reordering unit 40 comprises the output of apparatus ( see Fig.6b-c if 

20 X[n] is odd-length, or Fig.7d-e if sequence X[n] is even-length ). 

The block diagrams of apparatus implementing the subband 
reconstruction is shown in Fig.9. From the low-pass filter length and the high- 
pass filter length applied to the first 60 and second 62 input of the extension 
samples counter 64, it calculates the extension types for each channel and the 

25 number of extension samples k-low and k-high. The value of k-low coupled 
with the extension type are passed from the first output 74 to the first input 80 
of the low-band extension samples generator 86. At its second input 82 two P 
values of the starting/terminating points of symmetry is passed from the point 
of symmetry selector 78. The selector 78 retrieves those values from the starting 

30 and the terminating samples of the input sequence L[n], which is applied to its 
input. The P values are used only in the low-pass channel, because the points of 
symmetry in the high-pass channel are always equal to 0. The L[n] and H[n] 
samples are sent to the first 66 and second 68 input of the data reordering unit 
70. This unit performs inverse reordering ( see for example Fig.7f-g ). Then the 

35 reordered low-band samples via the first output of the unit 70 are passed to the 
first input 88 of the synthesis filter bank 108 and the third input 84 of the 
generator 86. The reordered high-band samples via the second output of the 



unit 70 are passed to the second input 90 of the synthesis filter bank 108 and the 
first input 100 of the generator 106. The extension samples are calculated 
separately in the low-band extension samples generator 86 and the high-band 
extension samples generator 106. Then the low-pass band extension and the 
high-pass band extension are sent to the third 110 and the forth 104 inputs of the 
synthesis filter bank 108 and together with samples from the first 88 and second 
90 inputs are used to perform signal reconstruction. The result of the 
reconstruction Y[n] ( see Fig.6f) from the output 112 of the synthesis filter bank 
108 is the output of apparatus. 

It will be appreciated that although only one particular embodiment of 
the invention has been described in detail, various modifications and 
improvements can be made by a person skilled in the art without departing 
from the scope of the present invention. 
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